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ABSTRACT 

We present the results from the discovery and study of ^ 3000 chromospherically active giant 
and subgiant stars toward the Galactic bulge. We find that these stars are predominantly RS CVn 
binaries with rotation periods between 10 and 100 days. We discover that the average rotational 
period of these stars decreases with their distance from the Galactic plane. We find that the 
primary stars in the RS CVn systems are predominantly first-ascent giants. Our research also 
suggests that, if these stars have spot cycles like the sun, then the cycle period must be longer than 
10 years on average. We confirm that the amplitude of the spot-induced modulations observed 
in the light curves of these objects is generally larger at minimum light than at maximum. 
Furthermore, we confirm that the amplitudes of the modulations due to stellar spots generally 
increases as the observed change in average brightness increases. We find no evidence for a 
relationship between a CA star's brightness and its rotational period. However, the average 
period does increase with colour for stars with periods < 30 days. 

Subject headings: Stars : Spots — stars : Variables: Other — Stars: Rotation 



1. Introduction 

Studies of chromospherically active stars (here- 
after CA stars) have been underway for many 
years. The two main types of spotted giant stars 
are single FK Com stars and binary RS CVn sys- 
tems (CAEs). FK Com stars are believed to be 
stars that have coalesced from binaries and gener- 
ally have shorter rotational periods than RS CVn 
binaries. Apart from the presence of strong Can 
H and K emission both types of CA stars exhibit 
photometric variability, soft X-ray and Ha emis- 
sion as well as UV and IR excess (Hall 1989). 
In many cases the secondary stars in RS CVn 
systems have orbital periods very close to syn- 
chronous with the primary star's rotation period. 
However, a significant number are far from syn- 
chronous (Strassmeier, Hall & Fekel 1993, Fekel & 
Henry 2005). To date only a couple of hundred 
of spectroscopically confirmed bright RS CVn bi- 
naries have been found (Strassmeier et al. 1988, 
1993). 
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The advent of large scale microlensing projects 
such as MACHO (Alcock 1995, 1997) and Ogle 
(Udalski 1994, 2000) greatly increased scientific 
knowledge about forms of stellar variability. Based 
on 4 years of photometry of Baade's Window, 
Olech (1996) found that a large fraction of the 
Miscellaneous Stars among the OGLE catalog of 
periodic variable stars were CA stars. However, 
this catalog was comprised of both ellipsoidal vari- 
ables and CA stars as a clear distinction between 
them was not found. Independently, in on-going 
studies into the nature of baryonic dark matter in 
the MACHO dataset we recently discovered the 
presence of a number of spotted variables toward 
the Galactic bulge and the Magellanic Clouds 
(Drake & Cook 2003). 

A number of detailed studies have been done 
on bright CA stars. These include determina- 
tions of the properties of both synchronous and 
asynchronous CABs derived from their orbital and 
photometric periods. Studies of their spots have 
progressed with surface images of the brightest 
stars. Spot locations have been derived from 
high resolution studies of individual spectral lines 
(Berdyugina, et al. 1999). In general the selec- 
tion of CA stars can be readily made by using 



the Can H and K omission strength from spectral 
surveys. However, the photometric properties of 
CA stars are also significantly distinct from other 
forms of stellar variability and they may be dis- 
covered from a few years of photometry. Clear 
distinction between RS CVn and FK Com types 
requires spectral information. In this paper wc are 
concerned with the general properties of CA stars 
on the Galactic bulge giant branch. We will use 
the large numbers of CA stars in this well defined 
region to examine the general properties of these 
stars. 

2. Observations 

The MACHO project repeatedly imaged ~ 67 
million stars in a total of 182 observation fields to- 
ward the Galactic bulge, the LMC and SMC, to 
detect the phenomenon of microlensing (Alcock, 
et al. 1997). These images were taken between 
1992 and 2000 on the Mount Stromlo and Sid- 
ing Springs Observatories' 1.3M Great Melbourne 
Telescope with the 8 x 2048^ pixel dual-colour 
wide-field Macho camera. In this analysis we will 
concentrate on the Galactic bulge where ~ 33000 
observations of ~ 33 million stars were taken in 
94 fields. Each observation field covers an area of 
43' X 43' on the sky. Most observations toward 
the Galactic bulge have exposure times of 150 sec- 
onds. The observations were taken with a dichroic 
to give simultaneous images in MACHO Bm and 
Rm bands. These passbands can be converted to 
the standard Kron-Cousins system V and R us- 
ing the photometric calibrations given in Alcock, 
et al. (1999). The conversion is made to Kron- 
Cousins V and R bands rather than B and R, 
as the transmission of the MACHO Bm band is 
more similar to a standard V band filters than a 
B band ones. Photometry was carried out using 
a fixed position PSF photometry package deriva- 
tive of the DoPhot package (Schechter, et al. 1993) 
called SoDOPHOT (Bennett, et at. 1993). 

The Galactic bulge was observed from March 
to October each year. The seeing of the data set 
was generally > 1" and measurements reach stars 
in each field of y ^ 21. The photometric sam- 
pling frequency varies greatly between the indi- 
vidual fields with the least sampled fields having 
only ^ 50 observations while those most sampled 
contain ~ 2000 observations. In most cases the ob- 



servations were evenly spaced during the observing 
season. Nevertheless, a small number of Galactic 
bulge fields have variations in cadence due to a 
changes in the observing strategy. 

3. Selection 

During our searches for Baryonic dark matter in 
the form of dark gas clouds we serendipitously dis- 
covered a number of CA stars toward the Galactic 
bulge (155) and the Magellanic clouds (12) (Drake 
& Cook 2003). Wc decided to undertake a specific 
search to select these variables from the MACHO 
database^. To achieve better selection we first ap- 
plied a number of cuts to the data to remove noisy 
or low signal-to-noise light curves. The initial step 
was to select light curves of stars exhibiting a sig- 
nificant amount of intrinsic variability. This quan- 
tity can be defined as: 

^^M^^^KM^, (1) 

V n — 1 n 

where by, is simply the weighted mean magnitude, 
and (Tft i are the mcasiired imccrtaintics associated 
with n individual magnitude measurements bi. We 
chose only light curves where a > 0.06 in both 
Bm and Rm passbands. As the stars were ob- 
served by the MACHO project simultaneously in 
Bm and Rm, we used the Welch & Stetson (1993) 
variability index: 

as a second discriminator for selecting the truly- 
variable objects. This statistic ensures that the 
variability observed in the two bands b and r is 
correlated. We required a value of / > 1.5, as 
used by Geha et al. (2003) for their QSO candidate 
selection from MACHO data. 

The candidates selected by the first two crite- 
ria include many different kinds of variable stars. 
In this analysis we are only interested in the CA 
stars. Therefore, the other variable star types are 
considered to he contaminants. In the analysis of 

•^data for variables stars from the MACHO project are avail- 
able at Web http;//wwwmacho.mcmaster.ca/. 
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Olech (1996) they noted the difhcuhy of separat- 
ing CA stars from ellipsoidal variables which occur 
in the same region of the CMD as CA stars. El- 
lipsoidal variables exhibit periodic sinusoidal vari- 
ability with a moderate amplitude similar to CA 
stars. However, ellipsoidal variables do not exhibit 
long timescale changes in brightness. In order to 
remove ellipsoidal variable stars from our candi- 
date CA stars we used the variogram function of 
Hughes, AUer, & AUer (1992) and Eyer & Genton 
(1999). 

The variogram can be determined by compar- 
ing pairs of measurements bi and bj taken at times 
ti and tj . We first determine n timescale bins over 
which differences in pairs of measurements are cal- 
culated by using: 



h{n)=nlog{0.1{to~tf)), 



(3) 



where to is the date of the first observation for the 
light curve and tf is that of the final observation. 
Within each time bin the median value, medp, of 
the difference between measurements squared is 
determined as: 



Varin) = \og{m.ed[{b, - bj)^]) 



(4) 



for bins h{n — 1) < \og(ti — tj) < h{n) and i ^ j. 
We only calculate the differences between mea- 
surements for observations separated by less than 
half the total time-span. If the entire time-span 
was used, there would be a strong bias in the me- 
dian value due to the decreasing number of pairs of 
measurements separated by increasing time differ- 
ences. We determine the slope of the variogram: 



5Var{n) 
5h{n) 



(5) 



for bins n = 3, 4, 5, as we are interested in bright- 
ness changes on long timescales. The pairs of mea- 
surement in these three bins are typically sepa- 
rated by times longer than 100 days. An increas- 
ing slope of the variogram selects objects that are 
variable on long timescales. For CA stars, changes 
in the mean magnitude occur over many years be- 
cause of the changes in the sizes and locations of 
star spots. By comparison ellipsoidal variables do 
not exhibit long term variations in baseline mag- 
nitude. For our selection we required that the can- 
didates have a modest increase in variability over 



time, Snar > 0.08, in each observation band. Eyer 

(2002) and Sumi et al. (2005a) used Syar > 0.1 and 
Svar > 0.09, respectively, for selecting QSO can- 
didates from the light curves of variable objects. 
Of the 33 million Bulge stars monitored the total 
number of candidates passing these three variabil- 
ity criteria is ^ 40000. This set of objects still 
contains most kinds of long period variable stars 
so further selections are required. 

Observations toward the Galactic bulge contain 
both disk stars and Bulge stars. It is difficult to 
determine the distances of the disk stars photo- 
metrically as they can occur anywhere along the 
line-of-sight. In contrast Bulge stars occur in a 
region < 1 kpc along the line-of-sight (Dwek et al. 
1995). In order to select a large set of CA stars 
at a uniform distance we decided to choose only 
CA stars near the Galactic bulge giant branch re- 
gion. In order to remove the disk dwarf stars it is 
necessary to correct the observations for the high 
reddening observed toward the Galactic bulge. We 
used the colour map of Popowski, Cook & Becker 

(2003) to correct for extinction in each Bulge field. 
As all the stars that we selected were variables we 
determined the median magnitude for each star 
and produce a CMD from the data. In Figure 1 
we present the colour-magnitude region selected to 
search for CA stars. This region contains ^ 11000 
variable stars. Stars to the left of the selection 
region are variable stars in the disk while those 
to the right are long period variables in the Bulge. 
The variable stars fainter than this box are a com- 
bination of faint Galactic disk stars and Bulge 
stars. The final step in our selection process is 
to search for periodic modulations due to spots 
among the variables in this region. 

4. Searching for Periods 

One useful property of CA stars is that their 
modulations can be used to determine their rota- 
tion periods. Modest changes in the spot induced 
period can occur because of the differential rota- 
tion of the stellar surface as well as the growth and 
decay of the spots (Fekel, Henry & Eaton 2002). 
However, the average periods derived over long 
timescales should reflect the rotation properties of 
these stars. 

For each of the stars in the selection region 
on Figure 1 we searched for periods between 2 
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and 150 days using the Lomb-Scarglc Discrete 
Fourier Transform periodogram method (Scargle 
1982) and Analysis of Variance (AoV) method 
(Schwarzcnbcrg-Czcrny 1989). Wc chose to only 
look for periods between 2 and 150 days as at 
shorter periods the daily observing schedule pro- 
duces strong aliases while additional aliases oc- 
cur at long periods due to the seasonal nature 
of the observing. Clearly this selection will miss 
CA stars with the shortest periods while at long 
periods we will show there are few such stars. 
In addition some of the CA stars are eclipsing 
while other are not. The eclipsing CA systems 
are likely to a different detection efficiency to non- 
eclipsing systems. We treated the Bm and Rm 
light curves both individually and in combination. 
We found the Lomb-Scargle periodogram method 
(hereafter LS method) to be more efficient at find- 
ing periods, which may be due to the sinusoidal 
nature of the variability in CA stars. The LS 
method also allowed us to determine the proba- 
bility that a period we found was due to a false 
alarm. For selection we required that there was 
< 1% chance of the period being due to random 
noise. However, we note the probabilities derived 
from the Scargle (1982) method were only taken as 
an approximation. More accurate confidence lev- 
els can be obtained via Monte-Carlo simulations 
(Schwarzenberg-Czerny 1996). 

Light curves that did not exhibit significant pe- 
riods were removed from consideration. Without a 
doubt there are a large number of CA stars among 
these as many of the objects removed exhibited 
long timescale changes in their brightnesses. The 
lack of a significant oscillation could simply mean 
that the spots are distributed evenly over the sur- 
face of a star during the observing period. A frac- 
tion of CA stars are also likely to be observed pole- 
on. These stars would display little modulation as 
the amount of observed spotted area would not 
change as the stars rotated. An excessive number 
of the candidates were found to have periods cor- 
responding to the moon's synodic and sidereal pe- 
riods. The periods of these objects were detected 
because of photometric problems caused by ob- 
serving fields near the moon. In most cases these 
objects were clearly not CA stars as they did not 
have a varying baseline and the periodic variations 
were not smooth. We removed these objects from 
our selection. In addition a large number of pe- 



riod aliases were discovered in field numbers 104, 
119, 301, 402. These were all observed many times 
per night during the last observing season. The 
change in sampling frequency led to a large num- 
ber of period aliases. Therefore, we excluded the 
data from the final season from the period search. 
As all of these fields have many observations this 
cut has little effect on our analysis. In Figure 2 
we present the V and R-band light curve of a CA 
star found in the Baade's Window region of the 
Galactic bulge. We did not make heliocentric cor- 
rections to the times of the observations as these 
are have little effect on the results. 

5. General Features of the CA stars 

We examined the light curves of all the can- 
didates by eye and removed all the light curves 
that appeared suspicious. During our selection we 
also checked accuracy of the periods by phasing 
the light curves of ^ 1000 of the light curves. To 
accomplish this we firstly fitted the light curves 
with a lower order cubic spline (polynomial and 
Fourier fits proved to be poorly constrained) and 
subtracted the fits. We then phased the residu- 
als with the periods found using the LS method 
(which was found provide accurate results). 

In Figure 3 we present the period distribution of 
the CA stars. The longest periods are around 110 
days. This fact suggests that even slowly rotating 
stars can have significant chromospheric activity. 
There are few CA stars with periods < 10 days 
where FK Com type variables are expected. This 
suggest that the bulk of the CA stars are RS CVn 
types. In addition, there appears to be some evi- 
dence of three peaks in the period distribution. 

After eliminating aliases due to the seasonal ob- 
servation pattern a number of the remaining ob- 
jects still exhibited multiple periods at the 99% 
level. We first investigated the most significant 
period then the subsequent periods. In most cases 
the additional periods were found to be due to 
aliases of a frequency peaks with very significant 
spectral power. However, it was clear that some 
of the periods were due to changes in the spots. 
In some cases the periods were twice, or half, the 
best frequency. Some were clearly due to a sec- 
ond group of spots appearing or disappearing at 
approximately 180 degrees from the original spot 
group. This gives rise to two modulations per ro- 
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tation. In most cases the longer period was clearly 
the rotational period when the light curve was 
phased. For some light curves we observed that 
the spot modulations shifted significantly with 
time. We believe this is related to the flip-flop 
effect where regions of activity disappear at one 
stellar longitude and appear at another (Bcrdyug- 
ina & Touminen 1998; Berdyugina et al. 1999). In 
Figure 4 we plot the V-band light curve for the 
object with Macho ID 109.19721.291. A smooth 
low-order spline was fit through this light curve 
and subtracted. The residuals have been phased 
with the 3.2 day period and these have been plot- 
ted in the lower two panels of this figure. The 
data points in the two lower panels have been sep- 
arated into the regions of different behaviour. The 
dashed lines in the upper panel mark the limits of 
these regions. In the lower left panel we plot the 
residuals of the points marked with dots in the 
upper panel. In the lower right panel we plot the 
residuals of the points marked by crosses in the up- 
per panel. Clearly there is a region of transition 
between the two frequencies at observation dates 
between 2400 and 2750 days since JD=2448623.5. 
Points in this region arc represented by triangles 
in the light curve but are not plotted in the lower 
panels. Clearly there are two real frequencies ob- 
served in this object's light curve. This change 
in the modulation frequency is likely to be due to 
the presence of a spot region located ~ 180° from a 
persistent group of spots (seen at phase 0.75 in the 
figure). As the main spot group recedes in size, or 
increases in temperature, the second group grows 
in size. This creates the modulation observed at a 
phase of ~ 0.25. The process appears to be cyclic 
as clear repartitions can be seen in the light curve. 

5.1. Eclipsing RS CVns 

As noted above all the light curves were ex- 
amined to remove uncertain candidate CA stars. 
During the examination we found a number of 
light curves that showed the signs of eclipsing bi- 
nary systems. We expect that any large group 
of binary stars will have some fraction that are 
eclipsing due to their random inclinations along 
our line-of-sight. The shorter the orbital period 
the more likely an eclipse will be visible. 

Measurements of eclipsing systems can be used 
to determine the orbital period and strongly con- 
strain the inclination of the system and the rela- 



tive sizes of the stars. If the reddening and dis- 
tances are known, as is the case for our bulge 
giants, accurate magnitudes of the stars can be 
determined and these can be compared with spot- 
less stars. As the spotted regions are eclipsed the 
times of spot eclipses can be used to determine 
the morphology of the spots. If high resolution 
spectral follow-up is taken, the masses, individual 
element abundances, and spectral types can be de- 
termined. This information can then be used to 
further constrain the radii of the stars and there- 
fore the absolute magnitudes and areas covered by 
spots. 

The discovery of eclipsing CA stars is not al- 
ways simple as in most cases the primary stars are 
giants of ~ IORq while the secondaries may be 
dwarf stars of ~ IRq. In such cases the change 
in total brightness due to an eclipse is small. In 
addition, the eclipse is superimposed on a vary- 
ing brightness level. Eclipses of synchronous RS 
CVn systems are not too difficult to discover as 
the eclipses occur at a constant location in the 
phased light curves. However, for asynchronous 
systems the eclipses occur at a different relative 
phase for each stellar rotation. Nevertheless, we 
did discover a number of eclipsing systems simply 
by examining the light curves. 

As an example of the eclipsing systems found, in 
Figure 5, we present the light curve, phased resid- 
ual curve and phased colour of eclipsing RS CVn 
system 108.19463.932 (MACHO ID). The RS CVn 
nature of this object is clear from the colour, mag- 
nitude, light curve modulation and change in mean 
fiux. We discovered this system has an eclips- 
ing period of 31.2 days. The flat bottom of the 
eclipse shows that secondary star is fully eclipsed 
for approximate 1.5 days each orbit. This eclipse 
is accompanied by a change in the colour showing 
that the secondary star is bluer. The depth of the 
eclipse is approximately 0.2 magnitudes indicating 
that it is fainter than the giant star. A small sec- 
ondary eclipse can also be seen in the phased light 
curve. As the phase separation of the eclipses is 
~ 0.5 the orbit is likely to be circular. The ex- 
act properties of the system depend on the actual 
masses of the components. However, the proper- 
ties derived from the light curves are consistent 
with a giant of radius 6.5R0 and a secondary with 
a radius of 2.OR0. The rotation period of the pri- 
mary appears to be ~ 31.6 days making this a 
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synchronous RS CVn system. 

5.2. Colours, Magnitudes and Locations 

We present colours and magnitudes of the CA 
stars discovered in Figure 6. In entire boxed region 
we found 2962 CA stars from among 11242 vari- 
ables. On this figure we also include an empirical 
estimate of the location of the Coronal Dividing 
Line (hereafter CDL, Haisch 1999) with the short- 
dashed line. The CDL marks the division between 
giants with hot coronac on the left and cool gi- 
ants with massive winds on the right (Henry et al. 
2000). Prom our sample, 56 CA stars lie to the 
right of the CDL where there arc 1638 candidates. 
There is little likelihood that these are not CA 
stars. Their locations in the CMD are probably 
due to them lying in regions where the amount 
of extinction is greater than that of the average 
values given by Popowski et al. (2003). 

In Table 1 we present the locations, V and V—R 
magnitudes, extinctions, changes in mean bright- 
ness AV and AR, and periods of the CA stars. We 
also mark the possible and definite eclipsing CA 
systems as well as the objects which show signs of 
the flip-flop phenomenon. 

The Galactic coordinates of the CA stars are 
given in Figure 7. The distribution of CA stars 
appears to follow the density of Bulge clump gi- 
ants, that is, increasing with decreasing Galactic 
latitude. However, at latitudes b > —2.5° there 
are few objects as the amount of extinction is very 
high. In addition, fields having b < —6° contain 
fewer observations than the other fields. There- 
fore, the sample of CA stars is likely to be less 
complete than lower latitudes. 

While investigating the period distribution of 
CA stars we discovered that there is a trend in the 
period with the Galactic latitude of the stars. In 
Figure 8 it appears that, on average, the CA stars 
closer to the Galactic plane have longer periods. 
To better illustrate this we have plotted the period 
density distribution for stars in the well sampled 
region between —5.4° < b < —2.1° in Figure 9. 
Here we binned the periods in 10 latitude bins by 
25 longitude bins. We then smoothed the distri- 
bution with a Gaussian kernel of 0.75 binwidth. 
This figure shows that, not only is there a clear 
shift in the peak of the period distribution with 
latitude, but the entire period distribution shifts. 



The slope of the contours make this point clear. 
The scale of the period change is approximately 4 
days per degree of latitude. It is not clear if the 
relation will hold at large galactic latitudes since 
this would shift the peak of the CA star period 
distribution to zero. 

The CA star catalog of Strassmeier et al. (1993) 
gives the locations of ~ 200 spectroscopically con- 
firmed binary systems. Almost all of these CA 
systems are within 200 pc of the Galactic plane 
suggesting that the objects in the catalog are pre- 
dominantly in the thin disk. At their observed 
scale heights these stars would have |6| < 2 at the 
distance of the Galactic center. In our experiment 
it is easier to discover CA stars with short periods 
near the Galactic plane than far from it because 
these fields were observed more often. Therefore, 
if a detection bias exists due to variations in the 
observing frequency it should lead to the detection 
of more short period CA stars near the disk. The 
observed period distribution shows the opposite 
eff'ect and is not due to such a bias. 

In Figure 10 we plot the cumulative period dis- 
tribution for the 558 CA stars in the region b > 
—3° and 650 stars b < —5°. We performed a non- 
parametric Kolmogorov-Smirnovtest on these two 
cumulative distributions. We find the greatest 
difference between the distributions {D = 0.224) 
occurs near a period of 20 days. The corre- 
sponding probability that these two period dis- 
tributions are drawn from the same population is 
P < 9.7 X 10~^^. This is very strong evidence that 
there is a change in the period distribution with 
Galactic latitude. 

In Figure 11 we present CMDs of the CA stars 
in the b > —3° and b < —5° regions. The stars 
in the high Galactic latitude region are fainter 
and bluer than those in the low latitude region. 
It is well known that the average age of Galac- 
tic disk stars increases with scale height. How- 
ever, it is not certain whether the Bulge formed 
rapidly or as the result of several accretion pro- 
cesses. Based on the analysis of RR Lyrae stars, 
Lee (1992) suggested that the oldest stars are in 
the central bulge region of the Galaxy. In this sit- 
uation one might expect that giant stars at low 
latitudes would be older and redder on average. 
The two CMDs in Figure 11 do appear to show 
this effect. Additionally, one expects that more 
evolved stars would rotate more slowly as they 
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proceed along the giant branch. The locations of 
the stars are consistent with this possibility. How- 
ever, Minniti et al. (1995) also found a relationship 
between the metal abundance of Bulge stars and 
their latitudes. This metalicity dependence com- 
plicates the interpretation or our result since the 
main-sequence turn-off magnitude is also depen- 
dent on stellar abundances. In addition, most of 
our CA stars are in binary systems. Such systems 
evolve differently than individual stars because the 
presence of a binary companion can delay changes 
in rotational velocity as a star proceeds along the 
giant branch. Contact binary star systems trans- 
fer mass from the primary star to its companion 
and this decreases the rotational velocity of the 
primary. However, due to their relatively long pe- 
riods we expect that most of the CA systems are in 
detached systems. For detached binaries the an- 
gular momentum loss of the system can cause the 
orbit to contract and system to spin-up (Maceroni 
1999). As the orbit spins- up the tidal interaction 
between the stars can cause an increase in the ro- 
tation rate of the primary. In this picture the more 
evolved RS CVn systems could have shorter peri- 
ods. We note that the increased rotation rate of 
these systems can lead to the increased magnetic 
activity seen in CA stars. Indeed, orbital periods 
changes have been observed for many RS CVn sys- 
tems (Hall 1989; Donati 1999; Lanza et al. 2001). 

5.3. Brightness and Period Variations 

In order to investigate the general properties of 
the CA stars we decided to look at relationships 
among their common parameters, such as peri- 
ods, dereddencd magnitudes, variations in mean 
brightness, and modulation amplitudes. 

The variation in the average magnitude of a CA 
star, AM, is a useful property since it reflects the 
change in the degree of spot coverage for a star. 
This value can be found by fitting the light curves 
of the CA stars. However, due to the large varia- 
tions in the number of observations, the periods, 
and the changes in average magnitude, as well as 
long gaps between observing seasons, we found 
it was not possible to obtain accurate values of 
AM from automated fits. Instead, we simply de- 
termined the difference between the 95th and the 
5th percentile of magnitude measurements. These 
measurement still reflect the variations of the light 
curve while rejecting the brightest and faintest 5% 



of points. Discarding the brightest and faintest 
points removes spurious observations and mini- 
mizes the effect of the sinusoidal modulations su- 
perimposed on the changing mean brightness. 

In Figure 12 we present the change in light in 
R-band, AR, versus that in V-band, AV. Clearly 
the ratio of these values is very well defined, al- 
though there are a small but significant number of 
outliers. Upon examination of the light curves it 
was clear that there were cases where bad data in 
the R-band observations lead to the outliers seen 
in the figure. Removing the outliers as shown on 
the figure we find that AR/AV = 0.900±0.068 for 
2703 stars, with uncertainty in the mean value of 
0.0013. This value is consistent with the result for 
RS CVn giants obtained by Henry et al. (2000). 
We note that the inclusion of the outliers does not 
change the result. We expect that uncertainties in 
the transformation from Macho Bm and Rm to V 
and R contribute uncertainty to this relation. 

In order to determine whether the amount of 
activity of the CA stars changes as they evolve 
along the giant branch, we decided to examine 
the relationship between each star's extinction cor- 
rected V-band magnitudes, Vq, and the change in 
its mean brightnesses, AV. To remove any bi- 
ases due to low detection sensitive for fainter stars 
we first excluded stars with apparent magnitudes 
V > 18. The left panel of Figure 13 shows there 
is a slight relationship such that fainter stars ex- 
hibit larger variations in their mean brightness. 
This suggests that faint stars have larger frac- 
tional spot coverage than brighter stars. There are 
two possible explanations of this result. Firstly, 
the faint stars have not evolved along the giant 
branch as far as the brighter stars, or secondly, 
these stars are simply fainter due to greater spot 
coverage. The observed average decrease in AV is 
only ^0.06 magnitudes over the range of magni- 
tudes displayed. To reduce the possible effect of 
an evolutionary change in magnitude along the gi- 
ant branch, we determined the ridge-line of the gi- 
ant branch and subtracted this from the observed 
magnitudes. The relationship between this resid- 
ual brightness, Vres, and the change in baseline 
magnitude is presented in the right panel of Fig- 
ure 13. After the ridge-line is subtracted there still 
remains some evidence for a relationship between 
the quantities. This result suggests that the trend 
in AV with Vq may simply be due to differences 
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in spot coverage, rather than due to evolutionary 
changes in the fraction of spot coverage. 

As a separate test we looked at how the proper- 
ties of CA stars vary with rotation period. Olech 
(1996) divided CA stars from ellipsoidal variables 
using "subjective judgment", and found that stars 
with periods < 10 days had an upper limit in their 
I-band magnitude. This limit was said to corre- 
spond to the limiting "break-up" period where a 
faster rotating single star would break apart. We 
found no evidence for a limiting V-band magni- 
tude for our CA stars. We note that long pe- 
riod variables (LPVs) show clear period magni- 
tude relationships (Cook et al. 1997; Wood 1999, 
2000). However, LPVs vary in flux due to pulsa- 
tion modes not due to surface features. We be- 
lieve that the lack of any such relationship in our 
analysis is evidence that our CA star selection is 
not contaminated with the ellipsoidal variables of 
Olech (1996). 

In order to further investigate the properties of 
the CA stars we analysed whether there was any 
relationship between the periods of the CA stars 
and their colours. To do this we binned the CA 
stars as a function of period and determined the 
mean colour of the stars in each period bin. For 
accuracy the bin size was adjusted so that there 
were 50 stars in each bin. In Figure 14 we present 
the colours of the CA stars versus their periods. 
The averages colours in each bin is shown as a tri- 
angle along with uncertainty error bars in these 
values. The average colour of CA stars appear to 
be significantly bluer at periods < 10 days than 
> 20 days. This may reflect the different nature 
of the two groups of CA stars (single FK Com 
stars and RS CVns). In the region between 10 
and ~ 30 days the CA stars become redder while 
at longer periods the average colour stays approx- 
imately constant. 

We also matched the CA stars with objects in 

the 2MASS all-sky point source catalog"*. We only 
included matches with measurements in J, H, and 
Kg bands and found 1077 matches with separa- 
tions < 1.5". The majority of the matched stars 
had J-K colours between 0.5 and 0.8 and J band 
magnitudes between 12.5 and 15. We corrected 
the near-IR values for extinction using the V-band 
values from the Popowski et al. (2003) 's map and 

^http://www.ipac.caltech.edu/2mass/releases/allsky/ 



the conversions given by Ricke & Lebofsky (1985). 
However, no firm relationship was found between 
the near-IR colours and the periods of the stars. 

5.4. Spots 

Rapidly rotating stars exhibit strong chromo- 
spheric activity. Therefore, one might expect the 
stars with the shortest periods to have more spots 
than stars with longer periods. We investigated 
this scenario by looking at how the change in mag- 
nitude (Ay) varies with the rotation period. As 
noted earlier IS.V is related to the fraction of a 
star covered by spots. We found that, if there is a 
relationship between these quantities, it is below 
the level of sensitivity obtained with this data. 

It has been found that many CA stars exhibit 
increasing modulations as their average bright- 
nesses decrease (Henry et al. 1995). However, 
there are also CA stars that exhibit very small 
changes in their mean brightness yet show largo 
changes in their spot- induced modulations (Henry, 
et al. 1995; Fekel, Henry & Eaton 2002; Fekel & 
Henry 2005) . The first effect is likely to be caused 
by increased spot coverage, since CA stars become 
fainter as spot coverage increases and the addi- 
tional spot area causes the modulation amplitude 
to increase. The second major cause is likely to 
be due to changes in the relative positions of the 
spots with time. The effect of differential rota- 
tion is well known for the Sun. For CA stars this 
effect could shift the relative positions of spots oc- 
curring at different stellar latitudes. As the spots 
line up on one side of a star they produce a change 
in brightness that increase with time, yet the av- 
erage brightness of the star remains unchanged. 
Other minor causes, such as varying spot tempera- 
tures and locations of the spots relative to our line- 
of-sight, can also cause changes in the observed 
brightness. 

In order to test how prevalent the effect of in- 
creasing spot coverage is on the modulation ampli- 
tude, we investigated the dispersion in each light 
curves at maximum light and at minimum light. 
Although the dispersion of the data points does 
not reproduce the actual modulation amplitudes, 
it does trace them. The dispersion value is less in- 
fluenced by outliers caused by bad data than the 
maximum and minimum data values. For sim- 
plicity we assumed that the data points followed a 
Gaussian distribution and we determined the asso- 



8 



ciatcd scatter about the mean measurement value. 
For each light curve we included all the data points 
within a two rotational periods of the maximum 
and minimum values. An exception was made for 
periods < 15 days where a 30 day period was used 
to insure a sufficient number of points to determine 
the dispersion. To further improve the accuracy of 
our result we removed CA stars with < 10 points 
in this time frame. 

In Figure 15 we present the difference between 
the dispersion in measurements taken near maxi- 
mum light, amax, and those taken near minimum 
light, (Jmim for V-band and R-band light curves. 
We found that 1317 of the 2186 (60%) light curves 
have ffmin — o-max > 0, where 1093 are expected 
if there is no relation between the dispersion and 
the stellar brightness. Therefore, as the modu- 
lation amplitudes directly follows the dispersion, 
the amplitude is greater at minimum light with 
4.8ij-level or with > 99.9% significance. The same 
result was also found for the R-band data. This 
shows there is little doubt that the largo modula- 
tions seen at minimum light are real. Nevertheless, 
in contrast, this figure also shows that there are a 
number of light curves with larger dispersions at 
maximum light than at minimum light. From our 
examination of the light curves it was clear that 
the maximum amplitude of modulation does not 
always occur at minimum light. There were light 
curves with small values of both amin ~ cTmax and 

To test how the changes in the relative posi- 
tions of the spots might cause large changes in 
the modulations we plotted the dispersion at mini- 
mum light versus the change in average brightness. 
The result is presented in Figure 16. 

There is a strong trend such that the greater 
the change in mean magnitude the larger the av- 
erage dispersion at minimum light. Few objects 
exhibit large dispersions without large changes in 
magnitude. This result suggests that the changes 
in the modulation amplitudes are primarily due to 
the changing spot sizes rather than the spots lin- 
ing up on the surface due to differential rotation. 

We attempted to model the evolution of stellar 
spots on a number of CA stars using the publicly 
available spot modeling code, such as SpotModelL 
of Ribfik, Olah & Strassmeier (2003). We found 
that this process was difficult for our data because 
the information about the spots is lost between 



the end of one observing season and the begin- 
ning of the next. The stars with long periods may 
be better modeled in this way as they will have 
changed less during the gap in observing. For the 
Magellanic cloud the MACHO project made con- 
tinuous observations throughout the year. We be- 
lieve these data would provide much more accurate 
models for the range of observed periods. 

5.5. Starspot Cycles 

One aspect of CA stars that has not been con- 
clusively answered is whether they have spot cy- 
cles similar to the Sun. A number of searches 
for these star-spot cycles have been carried out 
(Berdiugina & Tuominen 1998; Lanza et al. 1998, 
2001; Ribfik et al. 2003). In some cases cycles 
have been presented which correspond closely to 
the time-span of the data examined. Cases ex- 
hibiting only a single cycle might well be regarded 
with suspicion. Among the CA star light curves 
we observed, many objects passed through a pe- 
riod of maximum light more than once. Whether 
these truly correspond to a spot cycle is unclear 
to us. 

The amount of light observed from a CA star 
is dependent on the spot coverage, temperatures, 
and the locations of the spots. In general, we ex- 
pect that the minimum amount of spot coverage 
for a CA star will occur when it is brightest. The 
stellar longitudinal distribution of spots changes 
the amount of light observed over the period of 
star. Therefore, the mean light and spot coverage 
can be deduced by averaging the measurements 
over the star's rotational period. The movement of 
spots in stellar latitude relative to our linc-of-sight 
will also change the observed light curve, as will 
variations in spot temperatures. However, these 
effects appear less important than the changes in 
the areas of the spotted regions. 

In order to investigate the evidence for spot cy- 
cles and mitigate the secondary effects of spot tem- 
perature and stellar latitude, we decided to find 
the average of all the CA star light curves. Firstly, 
we shifted all the light curves in time so that they 
matched at their maximum light. Next, we sub- 
tracted the maximum light magnitudes from each, 
and combined the residual values. The time shift- 
ing process essentially gives us twice the time base- 
line of the observations since the time of maximum 
light is randomly distributed throughout the ob- 
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servation time. However, we did not include light 
curves where maximum light occurred within the 
first or last fifty days of observations, since we 
do not know whether wo have actually observed 
a maximum within the light curve in these cases. 

In Figure 17 we present binned residual CA 
light curve data points and the average light curve 
derived from these points for all the CA stars. As 
the time of maximum light was chosen to be within 
our observing season and there are large gaps be- 
tween observing seasons, the data points within 
the figure are unevenly distributed. In this figure 
each individual data point is the average of points 
within a 10 day-wide bin. For clarity, only 10 per- 
cent of the binned data points are plotted. The 
crosses show the median of the points at each time. 
These are representative of what one might expect 
the archetypal CA star might resemble in the ab- 
sence of spot-induced modulations. The median 
is presented in the figure as it better represents 
the points due to the long tail of faint values. We 
have only displayed the median values that derive 
from at least 15 data points and 10 difi'erent light 
curves. However, in most cases more than a thou- 
sand values were used to determine each median 
value. 

Figure 17 shows little evidence for the median 
magnitude rising to the value at maximum light. 
Although our CA star selection process required 
an increasing level of variability up to ~ 1000 days, 
our result indicates that average CA stars must 
have a spot cycle periods longer than 2000 days, if 
they indeed have cycles. There appears to be some 
symmetry in the shape of the data point distribu- 
tions on either side of maximum light and median 
values continue to decrease for ~ 1500 days on ei- 
ther side of maximum light. The combination of 
these two findings suggests that the average spot 
cycle should be longer than 3000 days. There is 
some evidence for an upturn in median brightness 
around 2000 days on either side of the maximum 
light. This may suggest an average spot cycle pe- 
riod of around 4000 days. However, there are far 
fewer points at the extremes of this figure, so this 
may simply be a numerical artifact. 

Both V and R band (not shown) light curves 
appear to be in good agreement. The V band 
data shows a larger amplitude as expected from 
Figure 12. It is possible that stars with difi'erent 
rotation rates will have different spot cycle times. 



To examine whether there were cycles visible for 
stars with similar rotation periods we investigated 
the shape of the median curve for subsets of rota- 
tion period. We found no significant evidence for 
differences from the combined data presented in 
Figure 17. 

6. Discussion 

In this analysis we have found approximately 
3000 CA stars. Without doubt our catalog does 
not include all of the CA stars in the observed 
fields. From simple geometrical considerations we 
expect to have found only a few of the CA stars 
which are pole-on to our linc-of-sight. For such 
stars a varying amount of spot coverage, or a 
change in spot latitude would result in a vary- 
ing baseline magnitude. However, these variations 
would not create significant modulations as the 
CA star rotated. Therefore, these stars will not 
exhibit the periodic modulations required for their 
discovery in this analysis. Additionally, as many 
CA stars may have many dozens of spots (Eaton, 
Henry & Fekel 1996), it is likely that some systems 
simply did not exhibit strong modulations during 
our observations because their spots were evenly 
distributed over the stellar surface. Furthermore, 
many CA stars in fields that were not observed 
frequently will have been missed. 

The CA giants we found in this analysis clearly 
occupy the same region of the CMD as clump gi- 
ant stars. Clump giants are used as standard can- 
dles for distances (Girardi & Salaris 2001) and as 
unbiased estimators of microlensing optical depth 
(Popowski et al. 2005, Sumi et al. 2005b). The 
presence of significant numbers of variable stars in 
this region of a CMD might bias the results derived 
from clump giants. We have discovered thousands 
of variables overlapping the clump region. How- 
ever, Popowski et al. (2005) found > 1.3 million 
giants observed in these fields, so the CA variables 
represent < 1% of stars in this region. In light 
of this, we expect that these stars will have little 
effect on clump distances or microlensing optical 
depths derived from clump stars. Furthermore, 
we find no evidence that the CA stars are concen- 
trated in the red clump region of the Bulge CMD. 
This strongly suggests that they are ascending the 
giant branch for the first time rather than clump 
giants. As most of our CA stars are located along 
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the giant branch, they are mainly within the Bulge 
rather than the foreground disk. Our discovery of 
a relationship between the rotation period and the 
Galactic latitude of CA stars could be due to vari- 
ations in the population with age and metalicity. 

One feature of CA stars is the presence of stel- 
lar flares. Flaring events can last up to 10 days 
but are relatively rare (Henry & Newsom 1996). 
Because our stars were observed simultaneously 
in two colours, we decided to search for outliers 
that might be due to stellar flares by looking for 
spikes occurring simultaneously in the two light 
curves. We found no evidence for stellar flares 
among the thousands of light curves. However, 
we expect that flares are much more likely to be 
prominent on dwarf type CA stars such as BY Dra 
types. 

Olech (1996) used low order Fourier fits in an 
attempt to separate ellipsoidal variables from CA 
stars. Olech (1996) 's method does not take into 
account the fact that the amplitude of the spot 
modulations changes with time or the fact that 
the baseline flux varies with time. We found their 
method to give very poor fits to our CA light 
curves. We believe this is the main reason Olech 
(1996) was unsuccessful at statistically separating 
ellipsoidal variables from CA stars. 

7. Future Work 

We discovered a number of eclipsing RS CVn 
systems and objects exhibiting the flip-flop effect 
by examining the light curves of the CA stars. The 
Fourier analysis of the light curves would provide 
a more complete set of such objects in our data 
set. 

In order to better understand CA stars it is 
necessary to study them in a variety of environ- 
ments. A study of CA stars in the LMC and the 
SMC would be enlightening, as the very limited 
analysis that was undertaken by Drake & Cook 
(2003) showed some evidence that there are fewer 
CA stars, per giant star, in the Magellanic clouds 
than the Bulge. The reddening toward the Mag- 
ellanic clouds is lower and better defined than the 
Bulge. In addition, the distances to the Magel- 
lanic clouds are well known. Furthermore, the 
Macho project continuously monitored stars in the 
Magellanic clouds during the term of the project. 
This coverage would eliminate period aliases seen 



in the Bulge due to the seasonal observing pattern. 
It would also allow the maximum and minimum 
light to be determined with greater accuracy than 
in this analysis. The analysis of Magellanic Cloud 
data could provide insight into the effects of age 
and metalicity variations on CA star periods. 

The spectral types of the CA stars discovered 
here could easily be determined using a few ob- 
servations with a wide field multi-object spectro- 
graph such as AAOmega^ on the 4 meter AAT 
telescope. This may provide insights into the na- 
ture of CA stars. The orbital periods of the sec- 
ondaries coidd bo determined for each of the CA 
stars presented. With this information the syn- 
chronously rotating RS CVn systems could be sep- 
arated from the asynchronous ones. In addition, it 
would be possible to test whether or not the rela- 
tionship we found between the period and Galac- 
tic latitude occurs only for synchronous or asyn- 
chronous systems. Synchronous CA systems may 
show trends in the rotational period with age as 
the strong tidal effects between the stars will dic- 
tate the observed rotational period. As the pri- 
mary evolves the redistribution of angular momen- 
tum is likely to change the orbital and rotational 
periods. Asynchronous CA star systems are more 
likely to have a rotation periods similar to individ- 
ual giant stars. However, this work would require 
high resolution spectroscopy taken on a number of 
occasions. 
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work. This paper utilizes public domain data ob- 
tained by the MACHO Project, jointly fimded by 
the US Departmcmt of Energy through the Univer- 
sity of California, Lawrence Livermore National 
Laboratory under contract No. W-7405-Eng-48, 
by the National Science Foundation through the 
Center for Particle Astrophysics of the University 
of California under cooperative agreement AST- 
8809616, and by the Mount Stromlo and Siding 
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Fig. Tlie light curves of one of the RS CVn 
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Fig. 4. — Tbii light (airv'cs of one of the CA stfirs 
exJiibiting (Jmngcs in period due to tht! (ycUc; nil- 
ture of two s|)ot regions at!p;u:.^te[i by 180" . See 
h!xt for further fJn^taiLs, 



Fig. 1. — Note the very poor figure quality on these figures is due to astro-ph's 1MB submission 
size limit. 
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Fig. 5. — Thp V-bfuxi light curve of an ptJijifiiug 
RS CVii st-iiT fouurl m tbe Galfictic ludge. The 
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Riihti^wting .1 spliiif; Bt tjo the light ciirvf?. The 
ookjiir trf the sj'steiui is; showo in the lo"wer light 
p^mel. 
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nmalj'SBd. The shdt-tt'Jshert line showg the aiipios- 
imate lc)c»ti(jn of the CDL. 
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Fig. S. — Tbe period dL^tribittioD of CA stJiis viei' 
siK their gaL^ctii; bititiicte. 



Fig. 2. — Note the very poor figure quality on these figures is due to astro-ph's 1MB submission 
size limit. 
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Fig. 9. — TLe smoot.bort deiEiitj' cfi!!tiibiitJou of 
CA star p«,iri(Mls vRisiis G.'ilactjc latitude in the 
weli sainpJetj regiou Ijc^tweeu — !5.4 < 6 < —2.1*. 
TLe CKiDtioiir JwtUs riispJiiycKl correspoKl to 10, 18 
and 2Ch stars. 





fig. 11. — Tltt^ (vloiiKi and inagiiitiidRi (rf CA 
stai? in Figure 10. Star» csEiuNticg Jong U'xni yivr'i- 
ability are shiTwn aa small dots wbik; CA stars are 
large (lota. In the left pani;l ive pr(5«!ii( tlie liwa- 
tioii tjf stius l> > —A" while iij. tlm riglit panel Tsp 
show CA stars h < -5°. 
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Fig. 10. — The cumulative praicHl (iistribiit ion tea 
CA shirs. Tbe snlid-liiu? re|a:(isc!iits CA st-ars ky- 
ftated at ^ > —3" whil^ the dasLed-line CA st,ars 
with b < -G°. 
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Fig. 12. — Tbe relatioi^liip best™™ the efiajige 
in V magnitiide versus R magmtudty liir nill CA 
stars . The objoct.s between the S ltt>rt-dashrai lines 
were used iu the calciil.atiou of the dis pension i£ 
the avTxage ratio of ampUt itctcs. 



Fig. 3. — Note the very poor figure quality on these figures is due to astro-ph's 1MB submission 
size limit. 
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Fig. 15, — TJie tiiangc in (iisp(3?!ioii hRt.wcssu miu- 
iiuiim ligjit iiud ijxiixi iDi uji light of tlie CA st^us 
VRigiiK tha dir,ivgif,'Mni at luiiiimiim light. V-b.-iu{{ 
d<ita. points are (lots aud H- bmd [X>iut5 iirp crosSfW- 



Fig. 14. — TIjb cciloiii of tbfi CA stjus veraiJH their 
pc!rio(ts. The triangles presont tint jii.Eriigfsi for 5(> 
objerts, wbikt thf. ra-ror bars give the uunertaiutiEW 
in each ai-eriige values. 



Fig. 4. — Note the very poor figure quality on these figures is due to astro-ph's 1MB submission 
size limit. 
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(AV. AR) 



Fig. 16. — The dispersioii in tlie brightuess of 
the CA st-axs fit minimiiui liglif ™rfnfs tbe rlmuge 
brigLtncsiis fiH' V jiriti R. Ugbt ciErves. V-Linct data 
points firo dots and R-baiKl points are crofsses. 




-30C0 -JtMW IDOO 2000 



AT Cdajfs) 

Pig. IT. — Tbe n^-ui brightxtesi of fdJ t be CA st-jvrs 
combined. Tbe cnKjses slioiv tbc uocxtiaD in efieb 
time iiiterMil. Only 10% of tbe points are sbown 
and eadi of ibeae represents tbe a'vijrage of over a 
10 day tiuu^spaii for e,acb ligLt curve. 



Fig. 5. — Note the very poor figure quality on these figures is due to astro-ph's 1MB submission 
size Hmit. 
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Table 1 
CA Star data. 



Macho ID 




RA 
h.m.s 




Dec 

o / // 




V 


V- 


-R 


A 


-V 


AV 


Ai? 


P 

(days) 


101.20649.254 


18: 


:04: 


:26, 


.39 


-27: 


:27: 


:44, 


.9 


16. 


,95 


1, 


,11 


2, 


.18 


0, 


,47 


0, 


,41 


67, 


.465 


101.20652.878 


18: 


:04: 


:31, 


.98 


-27: 


:19: 


:38, 


.0 


17. 


,66 


0, 


.86 


1, 


.88 


0, 


.29 


0. 


,25 


30, 


.754 


101.20653.340 


18: 


:04: 


:20, 


.86 


-27: 


:12: 


:51, 


.5 


17. 


,31 


1, 


.05 


1, 


.84 


0, 


.32 


0. 


,26 


26, 


.665 


101.20655.107 


18: 


:04: 


:25, 


.22 


-27: 


:04: 


:58, 


.4 


15. 


,99 


0, 


.89 


2, 


.04 


0, 


.22 


0, 


,20 


21, 


.803 


101.20657.810 


18: 


:04: 


:32, 


.36 


-26: 


:58: 


:44, 


.9 


18. 


,29 


1, 


.15 


2, 


.18 


0, 


.33 


0. 


,30 


49, 


.222 


101.20778.709 


18: 


:04: 


:49, 


.72 


-27: 


:35: 


:23, 


.0 


17. 


,50 


0, 


.97 


2, 


.04 


0, 


.23 


0. 


,18 


62, 


.050 


101.20779.1382 


18: 


:04: 


:45, 


.18 


-27: 


:30: 


:32, 


.8 


18. 


,77 


1, 


.02 


2, 


.07 


0, 


,27 


0, 


,20 


25, 


.306 


101.20780.380 


18: 


:04: 


:48, 


.23 


-27: 


:27: 


:32, 


.7 


17. 


,18 


0, 


.98 


2, 


.07 


0, 


.28 


0. 


,21 


59, 


.574 


101.20780.731 


18: 


:04: 


:47, 


.55 


-27: 


:24: 


:37, 


.6 


17. 


,70 


0, 


.91 


2 


.07 


0, 


.21 


0. 


,19 


33 


.427 


101.20780.840 


18: 


:04: 


:39, 


.84 


-27: 


:24: 


:31, 


.8 


18. 


,19 


1, 


,09 


2, 


.07 


0, 


,32 


0, 


,28 


62, 


.438 


101.20781.1471 


18: 


:04: 


:41, 


.27 


-27: 


:22: 


:37, 


.1 


18. 


,70 


0, 


.97 


1, 


.93 


0, 


.31 


0. 


,27 


14, 


.702 


101.20781.583 


18: 


:04: 


:36, 


.25 


-27: 


23: 


:08, 


.5 


17. 


,37 


1, 


.00 


1, 


.93 


0, 


.29 


0. 


,27 


51, 


.504 


101.20782.371 


18: 


:04: 


:45, 


.13 


-27: 


:18: 


:13, 


.7 


16. 


,93 


1, 


,02 


1, 


.79 


0, 


,35 


0, 


,30 


93 


.319 


101.20783.200 


18: 


:04: 


:47, 


.67 


-27: 


:12: 


:59, 


.7 


16. 


,59 


0, 


.90 


1, 


.85 


0, 


.70 


0. 


,64 


39, 


.191 


101.20783.854 


18: 


:04: 


:39. 


.29 


-27: 


:13: 


:54, 


.5 


18. 


,19 


0, 


.96 


1, 


.85 


0, 


.36 


0. 


,30 


52, 


.013 



Note. — Col. (1) gives the MACHO project object IDs. Cols. (2) & (3) supply the right 
ascensions and declinations at the J2000 epoch. Cols. (4) & (5) present the median magnitude 
and colour during observation period. Col. (6) provides the visual extinction derived from the 
Popowski et al. (2004) 's colour map. Cols. (7) & (8) show the maximum change in V band and 
R band magnitudes during the observing time. Col. (9) states the rotational period of the CA 
system. Table 1 will be published in its entirety in the electronic edition of the Astronomical 
Journal. The table is also available from author on request. 

^Possible eclipsing binary. 
^Definite eclipsing binary. 
'^Star which may exhibit the flip-flop effect. 
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